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Abstract 

We have estimated the contributions to the moments of polarized nucleon 
structure function gi{x, Q^) coming from the region of the very low x (10~^ < x). 
Our approach uses the nucleon structure function extrapolated to the region of 
low X by the means of the double ln'^{l/x) resummation. The Q^ evolution of gi 
was described by the unified evolution equations incorporating both the leading 
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order Altarelli-Parisi evolution at large and moderate x, and the double ln'^{l/x) 
resummation at small x. The moments were obtained by integrating out the 
extrapolated nucleon structure function in the region 10~^ < x < 1. 



1 Introduction 

The sum rules which are expected to be satisfied by the spin-dependent structure 
function gi of the nucleon play a very important role in the theory of the spin-dependent 
deep inelastic lepton scattering 0]. The sum rules involve (first) moments of the spin- 
dependent structure functions, and the moment integrals require knowledge of those 
structure functions in the entire region < x < 1 where, as usual, x denotes the 
Bjorken variable. Presently available experiments cover only the region of large and 
moderately small values of x (x > 5 • 10^^) for reasonably large values of the photon 
virtuality Q^ {Q^ > bGeV"^). Hence, a reliable theoretical estimate of the contributions 
to the moment integrals coming from the unmeasured small x region is important for 
the analysis of the sum rules. 

In this paper we propose an extrapolation of the spin-dependent parton distribu- 
tions and of the polarized nucleon structure functions into the low x region. The 
extrapolation is based on the double /n^(l/x) resummation |@,^. After integrating 
out the parton distributions or structure functions, one obtains the low x contribution 
to the corresponding moments. The integration interval extends from x ~ 10~^ to 
X = 1. It is assumed here that the small x behaviour of gi is controlled by the double 
ln'^{l/x) resummation. The full analysis of the double /n^(l/x) resummation effects 
was performed in detail in Ref. [0. The dominant contribution generating the double 
logarithmic terms is given by the ladder diagrams with the quark (antiquark) and gluon 
exchanges along the ladder. The very transparent way of resumming these terms is 
provided by the formalism of the unintegrated (spin-dependent) parton distributions 
which satisfy the corresponding integral equations. In |^ we extended this formalism 
so as to include the non-ladder bremsstrahlung terms by adding the suitable higher 
order corrections to the kernels of the corresponding integral equations. We also in- 
corporated the complete leading-order (LO) Altarelli-Parisi (AP) evolution within this 
scheme, thus obtaining the unified system of equations able to analyse simultaneously 
the parton distributions in the large and small x regions. In particular, this formalism 
allows us to extrapolate dynamically the spin- dependent structure functions from the 
region of large and moderately small values of x, where they are constrained by the 
presently available data to the very small x domain which can possibly be probed at 



the polarized HERA [|]. 

This paper is organized as follows. In section 2 we recall briefly Bjorken and EUis- 
Jaffe sum rules for nucleon structure functions. In section 3 the unified evolution 
equations for gii^XjQ"^) 0] which embody both complete LO AP evolution at large 
values of X and the full (ladder and non-ladder) double /n^(l/x) resummation at small 
X are discussed in the context of the partonic moment conservation. It is shown that 
in the non-singlet sector the first moments of both baryonic isovector g^^{x, Q"^) and 
baryonic octet g\ (x,Q^) are conserved, i.e. they are independent of Q^. It is also 
shown that there is no first moment conservation in the singlet sector. It should be 
recalled that the first moments of the non-singlet and octet structure functions acquire 
their Q^ dependence only as the result of the next-to-leading order (NLO) quantum 
chromodynamics (QCD) effects. Our formalism extends the LO AP formalism by 
including the small x resummation, yet it does not affect the conservation of the first 
moments of structure functions. 

In Section 4 our predictions for Bjorken and Ellis- Jaffe sum rules obtained after 
numerical integration of the respective nucleon components in the region extending 
from very low x (10~^ < x) are presented. First moments of nucleon structure functions 
are calculated and compared with experimental data fi^l^. In order to estimate the 



impact of the low x region on the sum rule integrals and moments, partial contributions 
from very low x region 10~^ < x < 10~^ are calculated explicitly. In our approach we 
use a simple semi-phenomenological parametrization of the non-perturbative part of 
the spin-dependent parton distributions. 

In Section 5 the summary of our results is given. 



2 Sum rules for gi[x^Q'^) 



The sum rules for polarized nucleon structure functions are derived from the space- 
time representation of scattering amplitudes T^f^[x) in terms of current commutators 
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ImT,k{x) = - (P, S I [3u3k]anUsym. I P' s) . (1) 

In the light cone limit x^ — *> 0, Xq — >■ 0, which corresponds to parton model kinematics, 
they reduce to: 



lim ImTik{x) = -eiki{p,s \ - 



J5/(0) + J^j|z(0) 



+ lj'iiO)\P,s). (2) 



The isospin symmetry determines the proton matrix element of the isovector current 
j|;(0), and results in the Bjorken sum rule for the non-singlet component of the nucleon 
structure functions gi'"'{x, Q^) which in LO approximation reads : 

^\^'ix,Q')=gA/6, (3) 



'0 

where gf'^{x,Q'^) = gi{x,Q'^) — g'^{x,Q'^), gl, g^ are proton and neutron structure 
functions respectively, and gA ~ 1.257 is the neutron /3-decay axial coupling constant. 
It should be stated clearly that the Bjorken sum rule acquires also corrections beyond 
the LO approximation. Since the formalism of the unified evolution equations we use 
henceforth includes only the LO AP evolution, we neglect the NLO correction terms 
both for the Bjorken and the Ellis- Jaffe sum rule. 

The Ellis- Jaffe sum rule for baryonic octet follows, if SU{3) flavour symmetry for 
octet /3-decays is assumed : 

^\!{x,Q') = {3F-D)/2A, (4) 



/o 
where : 

gl(x, Q^) = (Am + Ad- 2As)/24, (5) 

and F, D are octet /3-decay axial coupling constants [|1^] fulfilling the relation (3F — 
D)/2A ^ 0.024L Distributions Au, Ad, As denote quark components of the polarized 
nucleon. 

3 Moments oi gi{x,Q'^) and double ln^(l/x) resum- 
mation 

Low X behaviour of polarized nucleon structure function is influenced by double loga- 
rithmic Zn^(l/x) contributions, i. e. by those terms of the perturbative expansion, which 
correspond to the powers of /n^(l/x) at each order of the expansion p7|,p^. In what 
follows we will apply the double ln'^{l/x) resummation scheme based on the uninte- 
grated parton distributions [0,|T^,^. Conventional integrated spin-dependent parton 
distributions Api{x, Q^) {p = q, g) are related to the unintegrated parton distributions 
//(x', k"^) in the following way : 

AMx,Q') = Api'\x)+J^^ -^Mx' = x{l + —),e), (6) 

where Apl (x) is the nonperturbative part of the distribution, k'^ denotes the transverse 
momentum squared of the probed parton, W^ is the total energy in the center of mass 



W^ = Q^ (- — 1), and index / specifies the parton flavour. The parameter /cq is the 
infrared cut-off, which will be set equal to 1 GeV^. The nonperturbative part ApJ (x) 
can be viewed upon as originating from the integration over non-perturbative region 

Apr\x) = £'^Mx,e). (7) 

The nucleon structure function gi{x, Q^) is related in a standard way to the (integrated) 
parton distributions describing the parton content of the polarized nucleon : 



9'i{x,Q') = ^-^[9n^,Q')+9i'''ix,Q% (8) 

9ii^,Q') = ^-^[9U^,Q')+9i''''ix,Q% (9) 

where Nf denotes the number of active flavours {Nj = 3) and (e™) = — — ^^{ei)"^. For 

^f 1=1 
convenience we have introduced in (^, (^ the non-singlet and singlet combinations of 

the spin-dependent quark and antiquark distributions deflned for proton and neutron 

as: 



1=1 
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^^'^("H^, Q') = E 7^ - 1 ) (Agf^"^(^^ Q') + Agf "H^, Q')), (10) 



gf{x, Q') = E(A?7(^, Q') + Ag7(a;, Q')). (11) 

1=1 

In order to consider the Ellis- Jaffe sum rule we shall also analize the baryon octet 
structure function g^{x,Q^), deflned by equation (^. 

The full contribution to the double ln^{l/x) resummation comes from the ladder 
diagrams with quark and gluon exchanges along the ladder (cf. Fig. Q) and the non- 



ladder bremsstrahlung diagrams pl|,|22|. The latter ones are obtained from the ladder 
diagrams by adding to them soft bremsstrahlung gluons or soft quarks |1^, |T^, |21|, |22[ . 
They generate the infrared corrections to the ladder contribution. 

The relevant region of phase space generating the double ln^{l/x) resummation 
from ladder diagrams corresponds to ordered k'^/xn, where fc^ and Xn denote respec- 
tively the transverse momenta squared and longitudinal momentum fractions of the 
proton, carried by partons exchanged along the ladder p^,|2^. It is in contrast to the 



LO AP evolution alone, which corresponds to ordered transverse momenta. 

The structure of the corresponding integral equations describing unintegrated dis- 
tributions fNs{^', k^), fsi^', kP') and fg{x\ k^) for ladder diagram contribution read p[: 

hsix', e) = fSk^', e) + ^AP^M £ ^ /J'^ ^fNs K k'A , (12) 



'"mmm" 



'^mMm" 
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Figure 1: Ladder diagram generating the double logarithmic terms in the non-singlet 
component of the spin structure function. 



u2^ _ .(0). , ,,2^ , as r dz f^ '^ dk'^ 



fs{x',k') = f}^>{x\k') + 



2n 



/fc'2 



AP,,{0)fs{^,k'A+AP,,{0)fJ^,k 



12 



f,.,,>-rr'i^,,n^^i;/iC''-^ 



APg,{0)fs{-,k'- 



^Pggmg -,k 



with splitting functions APjj(O) = APij{z = 0) equal to : 



Apro) 



AP,,(0) AP,,(0) 
AP,,(0) AP,,(0) 



Nc ^ 



(13) 



(14) 



where as denotes the QCD coupling, which at the moment is treated as a fixed pa- 
rameter. The variables fc^(A;'^) denote the transverse momenta squared of the quarks 
(gluons), exchanged along the ladder. For the parton distributions in a hadron the in- 
homogeneous driving terms // (x', fc^) are entirely determined by the non-perturbative 
parts Apj- (x') of the spin-dependent parton distributions. 

Besides the ladder diagrams contributions, the double logarithmic resummation 
does also acquire corrections from the non-ladder bremsstrahlung contributions. It has 
been shown in Ref. p[ that these contributions can be included by adding the higher 
order terms to the kernels of integral equations (|12|) , ([T3|) . These terms can be obtained 
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Figure 2: The moment If of proton structure function gi{x, Q^) plotted as a function 
of Q^. Solid line denotes results obtained from the unified evolution including the 
full double logarithmic resummation /n^(l/x), dashed line shows pure AP evolution, 
dotted line corresponds to the non-perturbative input (overlaps with AP results). Ex- 
perimental data are denoted: E143 [0 with vertical bars, SMC |§] with stars, EMC 
with white squares. 



from the matrix: 
of the octet partia' 



Fg/^' 



[z)Go, where 



u^ 



z) denote the inverse Mellin transform 
wave matrix (divided by cj^ ), and the matrix Go reads: 



Gn 



2Nc ^ 

Nr 



(15) 



Following Ref. [^, we shall use the Born approximation for the octet matrix : 

T^Born' 

'z) = Air'^asMsln^z), 



UJ^ 



(16) 



where Mg is the splitting functions matrix in colour octet t-channel : 

Ms = 



1 

2Nc 


2 


Nc 


2Nc 



(17) 



In the region of large values of x the integral equations (|T2D, ([13|) describing pure 
double logarithmic resummation In (l/x), even completed by including non-ladder con- 
tributions, are inaccurate. In this region one should use the conventional AP equations 



25|^7| with the complete splitting functions APij(z) and not restrict oneself to the 
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Figure 3: The moment /" of neutron structure function (y'"(x, Q"^) plotted as a func- 
tion of Q^. Solid line denotes results obtained from the unified evolution including the 
full double logarithmic resummation ln^{l/x), dashed line shows pure AP evolution, 
dotted line corresponds to the non-perturbative input (overlaps with AP results). Ex- 
perimental data are denoted: E143 |10[ with vertical bars, SMC |jTl| with crosses, E154 
P] with stars, E142 [113 with white squares, HERMES [M with black squares. 



effect generated by their z —^ part. Following Refs. [@, 0, we do therefore extend 
equations ([T^), (|13]), and add to their right-hand-side (s) the contributions coming from 



the remaining parts of the splitting functions APij{z). We also allow coupling as to 
run, setting k'^ as the relevant scale. In this way we obtain unified system of equa- 
tions, which contain both the complete leading order AP evolution and the double 
logarithmic ln^{l/x) effects at low x. These equations are listed in Appendix A. 



3.1 Conservation of moments 



In order to get information about the moments of spin structure functions, we will follow 



the technique proposed in [0]. First we integrate the integrated parton distributions 
(B) over X : 



/ dxApi{x,Q'^ 
Jo 



'dxApr\x)+r-^ 



dxfi{x,e). (18) 



Let us denote Api{Q'^) = J^ dx Api{x,Q'^), and respectively fi{k'^) = Jq dx fi{x,k'^). 



Since structure functions Qi 



NS(S) 



X 



Q^) (p!OD, ([TI|) are linear combinations of (H), their 
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Figure 4: The ratio R^^ {Q^) for Bjorken sum rule integral plotted as a function of Q^. 
Solid line denotes results obtained from the unified evolution including the full double 
logarithmic resummation ln^{l/x), dashed line shows pure AP evolution, dotted hne 
corresponds to the non-perturbative input. 



moments may be obtained as : 



dk"^ 



9f{Q') 



_s,(o) , 
^1 + 



fcf) K yi -\- q2 



'JNsi^ )) 



k^ A;2(l + 



fc2 



fsik' 



where {i = NS, S) : 



_i,(0) 

^1 



dxg\{x), 
dxgi {x 



(19) 
(20) 



(21) 
(22) 



Furthermore, for non-singlet sector it was proven |T9| (see Appendix B) that the 
moments of fj^fsi^'i ^^) (0) vanish independently of the input gi (x) : 



Insk^ J 



0. 



(23) 



Therefore the unified equation for moments f^si^"^)^ obtained from equation (31) after 
integration over x, reduces to the integral equation with inhomogeneous term equal to 
0. Its solution then reads : 



JNsi^ 



0. 



(24) 



00 



0.014 
0.012 
0.01 
0.008 
0.006 
0.004 
0.002 






8 10 12 14 
Q^ [GeV] 



Figure 5: The ratio R^IQ"^) for Ellis- Jaffe sum rule integral plotted as a function of Q^. 
Solid line denotes results obtained from the unified evolution including the full double 
logarithmic resummation ln'^{l/x), dashed line shows pure AP evolution, dotted line 
corresponds to the non-perturbative input. 



For the singlet sector the situation gets more complicated. Although the quark 
singlet moment : 

(25) 

2\ 



f(o) 



fs'm = o 



vanishes again (see Appendix C), the moment of the input gluon distribution /^ (k 
takes a non-zero value which reads : 



p(0) 

g 



7T(k^) 



osik^ 

2tx 



2^f'(°) 



+ (^-f)A.f 



(26) 



Hence, the unified equations for quark singlet and gluon moments, obtained from cou- 
pled equations (0), (0) get a non-zero inhomogeneous term in the gluon sector, and 
the moment /^(A;^) is non- vanishing as well. 



The properties of partonic moments resulting from the unified equations 
(0) after transforming them into the moment space have clear implications for Bjorken 
and Ellis- Jaffe sum rules. Both Bjorken and Ellis- Jaffe sum rules (|^), (|^) concern 
evolution in the non-singlet sector. Hence, if one assumes input distributions g^ 
and g{ fulfilling the requirements (^, (^, it follows from relations (p!9|), (p^ ) that : 



9i 



■BJ^q2^ _ ^B.m 



9tiQ' 



9i 



= const, 
const 



(27) 
(28) 
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Figure 6: The ratio RP{Q^) for the first moment of proton structure function plotted 
as a function of Q^. Sohd hne denotes resuhs obtained from the unified evolution 
including the full double logarithmic resummation ln'^{l/x), dashed line shows pure 
AP evolution, dotted line corresponds to the non-perturbative input. 

are conserved throughout whole Q^ evolution. 

The non-vanishing unintegrated gluon input p6|) implies that there is no explicit 
conservation of ^f (Q^) ( and AglQ"^) ) during Q^ evolution. However, the conservation 
may be achieved implicitly by imposing the negative input gluon distribution Ap^^^ to 
fulfill the requirement : 



2^f'(°) 



+ ( 



11 



f)Apr 



0. 



(29) 



This is not a physical case but this shows that the moments are very sensitive to the 
gluon input which, in fact, has not an established phenomenological parametrization 
because of the lack of experimental data. Therefore, there is still possible to influence 
the Q^ evolution of ^i and henceforth, the evolution of T^'P'"'(Q'^) by manipulating 
the input gluon distribution. 
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Figure 7: The ratio R"'{Q'^) for the first moment of neutron structure function plotted 
as a function of Q^. Sohd hne denotes results obtained from the unified evolution 
including the full double logarithmic resummation ln'^{l/x), dashed line shows pure 
AP evolution, dotted line corresponds to the non-perturbative input. 

4 Numerical results for sum rules 



We solved the unified equations (|3T|) , (|32|) , (|33|) , assuming the following simple parametriza- 
tion of the input distributions: 



Ap: 



(0), 



x] 



mi 



X 



yVi 



(30) 



with r]u^ = rjd^ = 3, f]u = Vs = "^ and rjg = 5. The normalisation constants Ni were 
determined by imposing the Bjorken sum rule for Au^^ — Ad^^ and requiring that 
the first moments of all other distributions are the same as those determined from the 
recent QCD analysis [^. All distributions Api (x) behave as x" in the limit a; ^ 
that corresponds to the implicit assumption that the Regge poles which should control 
the small x behaviour of gi have their intercept equal to 0. 



It was checked that the parametrization (^) combined with equations (y), (|TOD, 
(|n|), (|3TD, (|32|), (p^) gives reasonable description of the recent SMC data on 5'f"^(x, Q^) 
and on gi{x,Q^) [0. After integrating out the respective parton distributions, we 
found that the moments ^f'^(Q^) and 'gf{Q^) are conserved during Q^ evolution with 
a good accuracy (not shown). The Bjorken sum rule is conserved explicitly, due to 
the choice of the input (y^ (Q^) ~ 0.2095). For the 'gf{Q'^) moment, the input was 
chosen to give ^i (Q^) ~ 0.016, which was implicitly in disagreement with hyperon 
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/3-decay data (H). 

We also investigated the Q^ evolution of the first moments of gi'"'{x, Q"^) and com- 
pared them with experimental data (see Figs. 0, ^. First moments of g\{x, Q^) agree 
well with the available experimental data both for moments obtained after performing 
LO AP evolution of non-perturbative input and for moments obtained after solving the 
unified evolution equation (cf. Fig. ^. On the contrary, our predictions for the first mo- 
ments of polarized neutron structure function gi{x, Q^) obtained from both the unified 
and the AP evolution are below the experimental data (cf. Fig. §). The discrepancy 
may be due to the fact that we consider only the LO AP evolution of the partonic 
moments. Also the Bjorken sum rule is considered at LO accuracy. Since the AP evo- 
lution dominates in the region of moderate and large x, applying it with the leading 
order (parton model) accuracy may be not sufficient to reproduce the experimental 
data. 

Moreover, we estimated the magnitude of contribution from the very low x re- 
gion (10^^ < X < 10^^) to the moments of polarized nucleon structure function. 
It was achieved by comparing the partial contributions with the total moments ob- 
tained after integrating g^'^ over x extending from 10~^ to 1. The calculated ratios, 

R = g\ (Q'^)/gUQ^) (z = BJ,8,p,n), where g'{^(Q^) = / dxgUx,Q^), are plot- 
ted in Figs. ^, §, ^, 0. For Bjorken integral (BJ) (|) the maximal | R^-^ \ is ~ 0.02, 
and for baryonic octet (8), | R |~ O.OL For proton (p) the maximal contribution of 
low X region to the first moment of the proton structure function is | i?^ |~ 0.02, and 
for neutron (n), | i?" |~ 0.08. 



5 Summary 

To sum up, we have estimated the contributions to the moments of polarized nucleon 
structure function gi{x, Q"^) coming from the region of the very low x (10~^ < x). Our 
approach used the nucleon structure function extrapolated to the region of low x by 
the means of the double ln'^{l/x) resummation which dominates in this region. The Q^ 
evolution of gi was described by the unified evolution equations, which incorporated 
both the LO AP evolution manifesting at large and moderate x, and double /n^(l/x) 
resummation dominating at small x. These moments were obtained by integrating out 
the extrapolated nucleon structure function in the region 10~^ < x < L 

Moments of proton structure function estimated for 2 < Q^ < 15 GeV^ were 
found in agreement with experimental data. The explicitly evaluated contribution 
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of low X region obtained via integrating out the gi{x, Q^) in the region of very low x 
(10~^ < a; < 10~^) was found to constitute about 2% of the total proton moment for 
2 < g2 < 15 GeV^. 

Moments of neutron structure function estimated for 2 < Q^ < 15 GeV^ were 
found to lie below the experimental data. This may be due to the fact that the unified 
evolution equations contain only leading order AP kernels. The contribution of the 
very low x region (10^^ < x < 10^^) was found to constitute around 8% of the total 
neutron moment for 2 < Q^ < 15 GeV^. 

Contributions from very low x region to Bjorken and Ellis- JafFe sum rules were of 
the order 1 % and 2 % respectively. 

This implies that the contribution of low x region enters the moments of polarized 
nucleon structure function at the level of 10 % at most. The contribution increases 
with increasing Q^. We expect that the improvements of the model needed to describe 
accurately the neutron data will possibly affect only the normalization of the corre- 
sponding integrals, and they will not change significantly the estimate of the relative 
contributions coming from the low x region. 



Appendix A 



The corresponding system of equations reads : 



(31) 



+ 



+ 



( Ladder ) 

asjk'') f^' d^A /■! dz {z + z^)/jV5(^, k'^) - 2zfNsix', k'^) 
Ikl fc'2 sL z 



27r 
27r 



fc' dfe'2 



fr2 



A;'2 



2 + -ln{l-x') 



l-z 

I lJ2\ 



fNs{x ,k 



( Altarelli — Parisi 




"m L 



^' dk'^ 



asjk^) /•! dz r^ /^ dk'^ 
27r Jx' z Jk'i 



fc'2 \^ 
(Non — ladder) 



fc2 fc'2 
qq 

•J2 



■fNs{-,k 



/2 



UJ^ 



k^^j2TT\ 



qq 



fNs[—,k 



/2 



fsix',e)=fl'\x',k') 



+ 



as{k'^) f^dz fi^y^dk'^A 



2n 



' z jq 



^'2 3^47'^' 
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+ 



+ 



2tt 
2tt 

2lT 

2tt 

asjk^) 

2tt 

as{k^) 

2tt 



1 dz /•^■'/^ 

' Z Aa ^ fc'2 



"'''fj^.'^' 



n 



( Ladder ) 

>^' dk'^ 4 /■! dz {Z + z2)/5(^, A;'2) _ 2z/5(x', fc'2) 



A;'2 3 



' z 



k' dk'^ 



fc2 k'"^ 



2 + -ln(l 
3 ^ 



I- z 

l2\ 



fs{x\k 






( Altarelli — Parisi 



1 dz 



w^ 



(^) 



27r2 



'^^ dA:'2 



gg '*=0 



A:'2 



/5 



1 dz /-^V^ dA:'2 



(Non — ladder) 



UJ^ 



^A;2^J27r2^ 



99 



U'2 



fA"-^k 



(32) 



/2 



/,(x',A;^)=/f(x',A;^) + 



+ 



+ 



+ 



+ 



asik"^) [^ dz f'^y^ dk 



■I z 



12 



i-2 



•:' Z Jkl k'"^ 

( Ladder 

^•' dfc'2 



/t'2 



'^^ Cifc'2 



2tt 

as{k'^) /■! dz /-^V^ dk 
2tt 

asjk') 
2tt 

asjk^) 

27r 

asik^) 

2tt 

as{k^) 

2-K 

asik^) 

2tt 



Wlf^iT-'" 



/2 



■12/, 



-,k'^ 

z 



/■I dz 4 fx' , 

Jx' 2: 3 \z 

fJ^-k'')-f,ix',k'') 



A;'2 



idz 

DZ 



'^^ dfc'2 



A;'2 
Altarelli 



11 iVi7 

- + 6ln(l 

2 3 ^ 



2/9 



jg(X , fc 



/2^ 



1 dz 

X' Z 



UJ" 



{< 



— Parisi ) 

27r2 



^ dz [^'I'dk'"^ 

x' Z 



fe2 A;'2 \^ 
(Non — ladder) 



59 

a;2 



'^^ dA:'2 



fc2 A;'2 



■/5 



z 



'^\G,\ 



A;2 



27r2/ -^nv'^' 

/ ( 



99 



z 



/2 



(33) 



In equations (|3lD, (p^), (^) we group separately: terms corresponding to the ladder 
diagram contributions to the double Iri^iXjx) resummation, contributions from the 
non-singular parts of the AP splitting functions, and finally contributions from the 
non- ladder bremsstrahlung diagrams. We label those three contributions as "ladder", 
"Altarelli - Parisi " and "non-ladder" respectively. 
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Inhomogeneous terms /j {x', k"^) 



NS, S, g), as stated above, may be expressed 



as : 



JNsi^ , k ) 



as{k-) 4 .1 dz (1 + z^)gr^'\^) - 2zgr^'\x') 



27r 3 Jx' z 



osik^ 

2n 



fPW,k') 



ff\x',k') 



2 + -ln{l-x') 
o 

2\ S{0)/^'x 



as{k^)A f^dz{l + z^)g', 



1- z 



-2zgf^'\x') 



27r 3 Jx' z 
as{k^^ « 



1-z 

(2 + ^Inil - x'))g',^'\x') 
Zn 3 



asik"^ 



"^dz 



271 



'-nJ ^{l-2z)Ag('\^), 

Jx' Z Z 



asik^)4: f^dz 



27r 



6 Jx' Z Z 

+ 6ln{l-x'))Ag'^°\x') 



2n 

2 3 



271 



6 



dz 



A^(o)(^)-2A^(o)(a;') 
l-z 



:i - '2z)Ag 



(0) 



,x 



(34) 



(35) 



Equations (|T]), (H), (||) together with (|3|), (H) and (D reduce to the LO AP 
evolution equations for nucleon structure function with starting (integrated) distribu- 
tions gi (x) {i = NS,S) and Ag^'^\x) after we set the upper integration hmit over 
k'"^ equal to k'^ in all terms in equations (0), (|32|), (p3D, neglect the higher order terms 
in the kernels, and set Q^ in place of W^ as the upper integration limit of the integral 
in eq. (§). 



Appendix B 



We prove that (^) holds. After integrating both sides of 
X = (0, 1) one arrives at : 



t(o) 



JNsik ) 



dx gi (x) 



+ 2 



as{k') 4 
27r 3 Jo 

«s(^^)_(o) , Ois{k'^ 



l-z 



over a; in the interval 



, l + Z^ /•! /l + Z^ 

dz h / dz 



l-z 



'9ns 



271 "''^ 27r 3 

Performing the integrals over z, one obtains : 

asjk^) 4 r a 
27r 3 



/ (ix ln(l — x) gi 
Jo 



NS,(0) 



X). 



f(o) 



JNS\k ) 



—2 / dx ln(l — x) gi 



NS,{0) 



X] 



+ 2 



as{k'^)^{o) , as{k 



27r 



^S^i? + ^27r 3 7 ^^ln(l-x)c/i 



<3_(0) 



x). 



(36) 



(37) 
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This implies : 



f(o) 



JNsi^ 



(3J 



Appendix C 



We prove that 
arrives at : 



7f(*^) 



(p^), (p^) hold. After integrating both sides of ( |35D over x = (0, 1) one 



S,{0), 



2tt 3 Jo 
:,2 



dx Qi (x) 



dz 

I — z 



1 (l + z^ 
dz 



\- z \- z 






5,(0), 



X) 



+ 



27r 



AT. 



1 dx 



_Ar,(0)( 



x^ 



Ap|,'^^(a;) / dz{x-1z). 



Using (|37| ) and performing integral over 2; in the last term, one obtains : 



In the gluon sector integration of (P^D over x yields : 



(39) 



(40) 



Ifi^) 



asik^) 4 /"i dx 



ZTT 3 JO X^ 



dz (2x — z) 



+ 
+ 



27r 



+ 1 - 2^1 - /" c/2 



"^^^^6 ['dxApf\x) \['dz(- 

Jo ^ VJo Vl — -2 

""^^^ ^6 C dx ln(l - x) Ap(°)(x). (41) 

JO ^ 



11 A^i.\^ 



Apf 



27r V 2 3 / "5 27r Jo 

Performing the integrals over z in (l4ll) , one obtains : 



7,0, . ^ a#^ 



^1 






(42) 
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